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Abstract. Let Cl C R 2 be a bounded domain with smooth bound- 
ary and a a smooth anisotropic conductivity on Cl. Starting from the 
Dirichlet-to-Neumann operator A„ on dCl, we give an explicit procedure 
to find a unique domain Cl, an isotropic conductivity a on fl and the 
boundary values of a quasiconformal diffeomorphism F : Cl — > f2 which 
transforms a into a. 



1. Introduction 

Let Q C R 2 be a bounded domain, and let a be a C 2 -anisotropic conduc- 
tivity denned over $7, i.e. a = (cr 13 ) is a positive definite symmetric matrix 
on Q in the C 2 class. The corresponding Dirichlet-to-Neumann map is the 
operator A a : C^dQ) -> C°(<9^) defined by 

(1-1) Kf = a ir \&a, 

where / G C 1 (90), ^ is the outer normal of and u is the C 1 (r2)-solution 
of the Dirichlet problem 

(1.2) V • (crVu) = on n, ulan = /. 



The equation (|1 .2[) represents the conservation of the electrical charge on 
SI if the voltage potential / is applied to d£l, and A a f is the current flux at 
the boundary. The following inverse problem arises from this construction: 
how much information about a can be detected from the knowledge of the 
mapping A CT ? 

It is not possible to determine a uniquely from A a . This was discovered 
by L. Tartar (see [7]). Indeed, let $ : f2 — > £1 be a diffeomorphism with 
&\dfi = Id, where Id is the identity map. Then we can define the push- 
forward of a as 

/'(Z>$)<r(Z>$)\ 
<!> ^' |det(S*)| J * ' 
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where D<5 is the matrix differential of <£, and one verifies that A$„ CT = A CT . 
In dimension two this is the only obstruction to unique identifiability of the 
conductivity. The anisotropic problem can be reduced to the isotropic one by 
using isothermal coordinates (Sylvester |1U|). and combining this technique 
with the result of Nachman for isotropic conductivities ([8]) we obtain the 
uniqueness result for anisotropic conductivities with two derivatives. The 
optimal regularity condition was later obtained by Astala-Lassas-Paivarinta, 
who proved the uniqueness for L°°-conductivities in [2]: for an anisotropic 
conductivity a £ L°°(Q) (Q C M 2 bounded simply connected domain) the 
Dirichlet-to-Neumann map determines the equivalence class of conductivities 
a' such that there exists a diffeomorphism : Q — > Q in the W ' 2 class with 
$\dn = Id and a' = <5*o\ 

The main purpose of this article is to clarify and show what one can 
explicitly reconstruct from a given Dirichlet-to-Neumann operator in the 
anisotropic case. From the results obtained in [10], [9], [8], [2], [6] we have 
deduced 

Theorem 1. Let Cl C M 2 be a bounded domain with C 1 boundary and let a 
be a C 2 -anisotropic conductivity on Q, isotropic in a neighbourhood of dQ. 
Suppose we know A$ : C 1 (dCl) — > C°(dCl). 

Then we can reconstruct a unique domain C I 2 (up to a biholomor- 
phism), an isotropic conductivity aonVt and the boundary values F\q^ of a 
quasiconformal C 1 -diffeomorphism F : Cl — > Q such that a = F*a. 

The new point in this statement is the existence of F : 0, — > 0, (and 
its explicit reconstruction at the boundary) without any assumption on the 
topology of f2. Early in [2j this result was proved for simply connected 
domains, a situation in which the question about deformations of complex 
structures of O does not make sense. 

Our main tool, as in |2], is the global solution F of a certain Beltrami 
equation equipped with an asymptotic condition, which takes our anisotropic 
conductivity a into an isotropic one, a, defined in general over a different 
domain Q = F(£l). With the help of F we then show the existence and 
uniqueness of a family of solutions ^){z,\) of the anisotropic conductivity 
equation, with special asymptotics at infinity, using also the existence of such 
type of functions in the isotropic case, that we call ip(w, A) (firstly introduced 
by Faddeev in [3]; see [9], [8] for the main properties). Then we show how 
one can reconstruct the boundary values of ip from the Dirichlet-to-Neumann 
operator A#, for any A, with a Fredholm-type integral equation, following 
the work of Gutarts ([6]). This is a generalization of R. Novikov's method 
for isotropic conductivities ([9]). We also show how to find the boundary 
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values of F from the knowledge of tP\qq (generalizing the result in [2]), and 
so we find F(dCl) = d£l (therefore also 0). 

After this, we explain how the knowledge of A#, iP\qq and F\ g ^ suffices 
to reconstruct the isotropic scattering amplitude 6(A). We give also another 
method: we define the anisotropic scattering amplitude b(\), and we show 
that it is equal to the isotropic one, proving that it is essentially a quasicon- 
formal invariant. This result was already included in [6]; here we give a new 
simpler proof. 

Thus with both methods, starting from ^P\qq we can reconstruct the 
isotropic scattering amplitude: this allows us to write the (9-equation which 
will permit us to find the isotropic conductivity a on 0, by the Novikov- 
Nachman reconstruction scheme ([9], |8j). 

Our scheme can be summarized in the following diagram 



A, 



Remark 1 . Although we cannot reconstruct a uniquely, for the applications 
it may be useful to find one representative of the equivalence class of a. To 
do this, using our theorem it suffices to find a diffeomorphism G : f2 — > f2 
with fixed boundary values (which are the boundary values of a quasiconfor- 
mal mapping, in our notation F\ d ^), and no other particular restriction: in 
this way (G~ 1 )*a will be a representative of a. If f2 is simply connected one 
can use the Ahlfors-Beurling extension theorem for quasi-symmetric home- 
omorphism of the circle ([IJ Thm. 2, p. 69]). 

Remark 2. We believe that an analogous result to our Theorem 1 is valid 
also on bordered Riemann surfaces. 




2. The Beltrami equation and Faddeev-type anisotropic 

solutions 

We identify M 2 with C by the map (x,y) i— > x + iy = z and we use the 
notation 

d z = ^(d x - id y ), dz = ^(d x + id y ) 

where d x = d/dx and d y = d/dy. We will also use the differential operators 
d, d such that df = d z fdz, df = d^fdz, with dz = dx + idy, dz = dx — idy. 
We also recall the identity d = d + d. 

We can suppose that a, already isotropic near 90, is the identity near d£l 
(see for the reduction to this case). Besides, we extend a to the whole 
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complex plane by putting a = I for z G C \ fi. Then, for the conductivity 
a = a 13 we define the following Beltrami coefficient 

, , -a ll {z)+a 22 {z)-2ia l2 {z) 



d- ll {z)+a 22 {z) + 2^/d^) 

which satisfies |^i(z)| < k < 1 and is compactly supported in £1. We now 
recall the existence of a diffeomorphism that transforms a into an isotropic 
conductivity. 

Lemma 2. (Sylvester |10| ) There is a quasiconformal C 1 -diffeomorphism 
F : C — > C suc/i that 

F(z) = z + O ( - ) as |z| — > oo, 



and /or which 

(F*a){z) = a(z)I := (det(o-)) 1 / 2 o F~ 1 {z)I. 

Thanks to results by Ahlfors and Vekua ([1], .F is obtained as the 

solution of the Beltrami equation c\F = [Aid z F, so F is holomorphic in C\0. 

Proposition 3. There exist unique Faddeev-type solutions of the anisotropic 
conductivity equation, i.e. functions ip(z,X) such that 

(2.1) V • (o-(VV-)) = 

for all z G C, AG C, and tp(z, A) = e Xz (l + 0(-)) when z — > oo. 

Proof. We define = i^(f2) and g = ^rpr- It is known that if u is a solution 
of V • (ffVu) = in fi, then u = (T 1 / 2 ^ is a solution of 

(2.2) -An + <?tt = 

in $7. From [3] and [8j Thm. 1.1], we have that for every A G C there is a 
unique solution tp(w,X) of (|2.2p with the asymptotic behaviour ip(w,\) = 
e Xw (l+0(^)) when w — > oo. So we directly have that ip{w, A) := a~ 1 ^ 2 ^jj(w, A) 
is a solution of V • (crVi/O = with the same asymptotic (because a = 1 
outside Q). 

Now let i/j(z, A) be a Faddeev anisotropic solution. If we consider ijj'(w, A) = 
ifj{F~ l (w), A), we have that V-(<rVV' / ) = from the construction of a. Using 
the properties of F and tp, we get, for w — >• oo, 



f(^,A) = HF-Hw), A) = e A ^» ^i + f_I 
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showing that ifj'(w,X) satisfies the same asymptotic of ip(w,X). From the 
uniqueness of ip(w, A) we obtain 

(2.3) i>(z,X) = i;(F(z),X), 

which proves both existence and uniqueness. □ 

From the equality (|2.3p we can also derive a useful formula to calculate 
F\qq. In fact, results in [5j also indicate how the family of Faddeev solutions 
behaves with respect to A. We have indeed \e~ w ip(w, A) — 1| — > as |A| — > oo 
for every fixed toeC. If we take w G C \ O the same limit is valid also for 
ip(w, A); combining this with (j2.3|) we deduce the following formula. 

Proposition 4. For all z G C \ fi (in particular for z G dCl) we have 

lim logWM)) = lim hgW(fW,A)) = F(z) 

|A|— >oo A A|— >oo A 

3. An integral equation for ip\ d ^ 

Following the approach of [6], we show that, as in the isotropic case, we 
can find iP\qq through a Fredholm-type integral equation. 

The main idea is to decompose the differential operator —V ■ cxV as 
—A + M, where M is a compactly supported operator. So we can char- 
acterize ifj(z,X) as the solution of the following integral equation: 

ip(z, A) = e zX — - / G(z - w, X)Mip(w, X)dw A dw, 
2 Jn 

where 



^. >N ie I e 2 aw A aw .. , 

G(z, A = / , z G C, A G C 

is the Faddeev-Green function for the Laplacian. 

Proposition 5 ([6]). For every A G C the boundary value of tp satisfies 

(3.1) 4>(z,\)\ dCl = e zX - / G(z-w;,A)(A^ - Aq)^^^)^, 

where Aq is the Dirichlet-to-Neumann operator of the standard Laplacian ( or 
for the case of constant conductivity). 

This follows from the identity 



(3.2) / u (Aa-A )u= / u Mu, 

Jan Jn 

where u ,u G W 1 ' 2 ^), V • (aVu) = 0, Au = in 0. 

The fact that the integral equation (|3.ip is of Fredholm type in the Sobolev 
space W s ' 2 (dCl) is the content of [6l Lemma 2.5], and the invertibility is based 
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on the absence of exceptional point for the equation (|2.ip (which is implied 
by the same result in the isotropic case) . 

4. Reconstruction of the scattering amplitude 

Following [5], we define the non-physical scattering amplitude for the 
isotropic inverse problem as 

(4.1) 6(A) = / e- Jw q(w)4>(w,\)dw. 

From [9j we have 

6(A) = f e- Jw (A q -A )i>(w,\)dw, 
Jan 

where A q is the Dirichlet-to-Neumann operator of the Schrodinger equation 

{22}. 

Since a is the identity near dQ, equation A q = a~ 1 / 2 (A r7 + \^) ~~ 1 ^ 2 
reads A q = A a , and 4>\dn = t/j\dQ, so 



(4.2) b(X)= e- Xw (A a - A )iP(w,X)dw. 

JdCl 

Thus, for the reconstruction of 6, it is sufficient to determine A a and ip\an- 
By (|2.3p we already know ip\dQ] for the determination of A a , by arguments 
of [2j, we obtain the identity 



(4.3) / uA a v = / uA a v 

Jdh Jan 

which holds for any u, v G C 1 (c?Sl) and u, v G C 1 (dQ) such that u = uoF and 
v = v o F (this follows directly from the properties of F and the symmetry of 
the two Dirichlet-to-Neumann operators). So we find A a from A# and F\qq- 

4.1. Complementary results. Here we give here another method to find 
6(A). Inspired by [B], we define the anisotropic scattering amplitude as 

(4.4) 6(A) = [e- J *M4>(z,\)dz 

and we have the following result. 
Proposition 6. 6(A) = 6(A) 
Proof. From identity (|3.2p we find 



6(A) = / e- Az (A a -A )i;(z,\)dz. 
I an 
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Using the relation f|4. 3|) (that is trivially true also for Ao), if we call G 
we arrive at 

6(A) = / e" AGM (A ff - Ao)#(G(w), X)dw 
Jan 

(4.5) = / e- XG( - w \A a - A )iP(w,\)dw 

also by (|2.3p . Now we only have to prove that (|4.5p is the same of (|4.2p ; we 
will use the following lemma. 

Lemma 7. For every tfi € C 1 (dQ), tp E C 1 (il) solution of V • (aVip) = in 
£1, we have 

(4.6) I <t>(K - A )V = 2i [ 0(00 - 500 ), 

where AtpQ = in Q and ipo\dU = ^IdQ- 

For now we assume this lemma, and we show how it takes us to the 
conclusion. 

Let Or = { z £ C : \z\ < R } the disk of radius R, and let R be 
sufficiently large to have SI C We extend 0o to the whole plane and 

we apply Stokes' theorem to f2/j \ f2. We obtain, for every quasiconformal 
homeomorphism E : C — > C, holomorphic in C \ 0, 

/ e -^H(^_^ ) = / e -AB(»)(^_^ fl ) 

Jan ./an^ 

+ f a( e -^)(5v-^o)) 



but the last term vanishes, because e ^ E ( w ) is anti-holomorphic and 300 = 
in C \ O. So the identity 

[ e- XE(w) (d^-diJo) = [ e- XE(w) (dij-d4jo) 

Jdn Jon R 

is true for R ^> 0, E(w) = G(w) and E(w) = w. As we have G(w) = 
w + O(t^t) for w — > oo, using the lemma we deduce 
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b(X) = 2i / e- XG( - w \di> - dipo) 
Jan 



lim 2t / e- XG ( w \diJ> - d^ ) 



an f 



lim 2i [ e~ Xw (di(}-dipo) 
Jan* 



= 2i f e-^(d^ - d^ ) = 6(A) □ 
Jan 

Proof of Lemma 7. Let a € C 1 (Q) such that a\gQ = 0, and w = x + iy. 
As we have already remarked that A CT = A q , it suffices to prove the lemma 
for a solution i/j of (|2,2p . From the definition of the Dirichlet-to-Neumann 
operator and from Stokes' theorem, one has 

/ (f>(A q - A )V> = / (Va • V(V> - Vo) + qaip)dxdy, 
Jan Jn 

and by Stokes' theorem and by the identity A = 4^^= 

2i / 4>(dip — dipo) = 2i da A (dtp — dtpo) + / Qa>i> dxdy. 
Jan Jn Jn 

Writing in coordinates we get 

da A ONj — ifro) = — Va • V(ifi — if>o)dxdy + -da A d(ip — ipo). 
2i 2 

Again by Stokes' thorem we have 

da A d(if) — ipo) — — I (tp — ipo)da = 
n Jan 

because tp\ga = V'olan- The proof follows. □ 

5. The 9-equation and the reconstruction of a 

Here we follow the steps of |9j to reconstruct isotropic conductivities. 
The function fi(w,X) = tp(w, X)e~ Xw satisfies the following 9-equation with 
respect to A 

(5.1) = b Me^]^\). 

d\ 4irX 

This is equivalent to the integral equation: 

(5.2) u(w, A) = 1 + — ^ [ b ^ _, e J ' w - x ' w fi(w,X')dX' A dJ' 

8tt 2 z 7 c (A' - A)A 
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because \x — >• 1 when w — > oo. By results of [8], equation (|5,2p is solvable, 
and one can find a(w) from the integral formula 

(5.3) ^^(w) = 1 + — ^- / ^e^- Aw >(«;,AWA A dA, \/w € C 

8tHz 7 c |A| 2 

or from the more stable formula 

(5.4) = — — , Vw € C, VA € C. 
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